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Global  IFS  seem  to  be  suited  best  for  compressed  encoding  of  natural  objects  which  are  in  most  cases 
self  affine  even  if  not  always  exactly.  Since  affine  transformations  - the  IFS-Codes  - resp.  the  union  of 
all  their  orbits  generate  an  object  (an  IFS- Attractor),  the  detection  of  a non  minimal  set  of  these  orbits 
solves  the  inverse  IFS-Problem  by  calculating  a superset  of  IFS-Codes,  which  has  to  be  minimized. 

Here  a method  is  presented  how  these  orbits  (in  particular  those  on  the  object  boundary)  can  be 
calculated.  Therefore  a generalized  convex  hull  - the  rr-Hull  - is  defined.  This  fractal  hull  is  bounded 
by  log  spirals,  that  curves  formed  by  the  orbits.  It  is  shown  that  log  spirals  can  be  represented  by  a 
continuos  function  of  powers  of  affine  maps  and  that  by  using  this  "spiral  equivalent"  the  generating 
transformations  of  the  orbits  by  which  an  IFS-object  is  bound,  can  be  calculated  in  the  x/y-plane. 
Further  more  this  representation  can  be  used  for  the  classification  of  the  detected  orbits,  necessary  to 
calculate  the  IFS-Codes  of  a minimal  IFS  from  their  generating  transformations,  subsequently. 


1 Introduction 

1. 1 Problem  Specification  - Objects  to  be  encoded 

The  idea  is  to  globally  encode1  objects  of  an  image  appearing  in  the  nature  and  therefore 
often  self  similar  to  a high  degree  but  difficult  to  encode  or  to  compress  by  a 
mathematical  formula.  This  global  IFS-Approach 2 - representing  an  object  by  the  union 
of  affine  contractive  transformed  copies  (subobjects  generated  by  the  IFS-Codes)  of  the 
object  itself  - seems  to  be  used  best  to  solve  the  inverse  (global)  IFS-Problem. 

Thus  the  Collage  Theorem  which  can  be  used  to  generate  whole  objects  of  an  image 
is  given:  A set  {cn(:  R2  R2 1 affine,  contraction)  of  complete  Euclidean  space  (R2,d) 
is  an  Iterated  Function  Systems  (IFS)  with  conductivity  factor  s,  for  i e N,  the  distance 
d,  0 < s < 1 and  s-h(B,  C)  > h(W(B),  W(C)),  if  for  the  Hutchinson  operator  W°  holds 

W( A)  = U U ®i  ( x)  and  A = W( A). 

cqew  I6A 

The  compact  set  A is  called  an  attractor  of  W. 

The  Collage  Theorem  (CT)  for  attractors  close  to  given  compact  sets: 

h(L,  W(L))  < e =>  h(L,  A)  < e/(i.s) 
where  h is  called  the  Hausdorff  Distance 


a W:  rSfR2)  ->  WR2)  with  W(A)  = (W(a)  | ae  A),  Ae  ’TAfC)  is  a contractive  map  of  (^R2),h)  the 
space  of  nonempty  compact  subsets  of  R2  where  h(C,B)  = max{d(C,B),  d(C,A)}  and  C,  Be'JWR2). 
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L c R2  is  a given  compact  set  of  (R2,d)  and  s > 0.  For  e = 0 (L  = A)  the  CT  for  exactly 
self  affine  attractors  is  obtained.  In  principle  the  CT  can  be  applied  to  any  object  if 
sufficient  many  IFS-Codes  ©;  are  used.  However,  to  compress  an  object  the  number  of 
parts  (L)  (the  subobjects)  have  to  be  small  for  a small  s. 

For  the  discrete  pixel  space  the  discrete  CP  is  defined  as  a set  {©m:  P2  ->P2 1 wm 
affine,  contraction}  of  the  pixel  space  (P\d)  - P = (0,...,R-1 } c N - is  a Discrete  Iterated 
Function  Systems  (DIFS)  with  conductivity  factor  s,  for  m e P,  R2  the  number  of  image 
screen  pixels,  the  block  distance  d,  0 < s < 1 and  s-h(B.C)  > h(W(B),  W(C)),  if  for  the 
Hutchinson  operator  W:  *?<(P2)  -»  ^P2)  holds 

W(A)  = (J  (JmJ*)  and  A = W(A). 

meP  xeA 

The  set  A is  called  an  discrete  attractor  of  W.  The  discrete  attractor  Ae^P2)  is  not 
unique,  but  there  always  exist  a unique  maximal  attractor  Aj^e^fP2)  with  A= W(A)  => 
AcA,,,^  which  contains  all  discrete  attractors. 

The  Chaos  Game  generates  an  attractor  A by  the  iterated  random  application  of  the 
different  com  to  one  image  point  i.e.,  by  the  union  of  orbits  formed  by  all  combinations  of 
®m  - the  basic  principle  of  the  presented  solution. 

In  contrast  to  global  IFS,  an  encoding  based  on  Partitioned  IFS  (PIFS)4  divides  the 
image  into  larger  and  smaller  segments  (polygons).  The  smaller  target  segments  are 
interpreted  as  affine  contractive  copies  of  the  larger  ones  and  then  the  CT  is  applied. 
Since  only  IFS-Properties  of  affine  relationships  constrained  to  the  segment  form  can  be 
exploited,  all  of  the  many  image  parts  are  only  approximated  by  an  IFS-encoding. 

1.2  The  Solution  - How  to  encode  global  IFS-Objects 
Calculating  boundaries  of  an  object 

In  general  the  boundary  of  an  object  contains  the  most  information  how  an  object  is  built. 
Points  inside  an  object  normally  cannot  be  related,  because  they  are  difficult  to 
distinguish.  The  common  method  to  determine  the  boundary  of  an  object  - especially  if  it 
is  defined  in  a discrete  (pixel)  space  - is  to  form  the  hull  of  an  object. 

We  first  can  try  it  by  the  convex  hull  #A),  which  is  defined  as  "Intersection  of  half 
planes  containing  a point  set  (the  attractor)  A".  It  soon  appears  that  a hull  of  such  a kind 
is  not  usable  for  IFS-Attractors,  because  only  the  boundary  of  convex  objects  - what  an 
IFS  certainly  is  not  - can  be  determined. 

Second  we  can  try  a more  general  hull  - the  a-Hull5  >?a  - defined  as  follows  : 

The  a-Hull  ^a(A)  (see  Fig.  1)  is  defined  as  the  intersection  of  all  generalized  disks  B(c,  r) 
resp.  B(c,  -r)  of  non  negative  radius  r = !/a  and  the  centre  c containing  the  attractor  A, 
where  B(c,  -r)  is  the  complement  of  the  disk  B(c,r).6 

The  boundary  of  an  a-Hull  is  formed  by  a-Extremes  (see  Fig.  1): 

A point  s e A is  termed  a-Extreme  (comer  of  *?a(A))  in  A,  if  there  exists  a B(c,  Vo) 
such  that  s lies  on  its  boundary  and  contains  all  other  points  of  A.6 

Now  it  is  possible  to  analyse  concave  objects  and  even  objects  concave  with  different 
curvature  using  different  radii  of  the  generalized  disks.  But  we  need  a relationship  of 
boundary  points  (extremes)  which  provides  more  information  about  the  object.  Using  the 
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a-Hull  we  get  only  a relationship  of  points  laying  on  the  same  edge  of  a disk  of  different 
radius  forming  the  a-Hull.  But  the  boundary  points  of  an  IFS-Attractor  do  not  lay  on 
circles.  They  form  orbits*  in  the  discrete  space  which  are  given  as  subset  of  W+  c for  an 
IFS  W (as  it  is  easily  to  observe  in  the  application  of  the  Chaos  Game).  These  orbits  lie 
on  logarithmic  spirals1,  what  will  be  shown  in  this  paper.  So  the  best  suited  hull  to 
determine  an  IFS-Attractor  seems  to  be  the  a-Hull  which  is  bounded  by  log  spirals  and 
will  be  detailly  defined  and  described  in  the  main  part  of  this  paper. 

Calculating  a minimal  IFS 

Once  found  the  orbits  on  the  boundaries  of  different  a-Hulls  they  have  to  be  classified  to 
calculate  the  IFS-Codes: 

Classes  C;  = {orbits  which  can  be  affinely  mapped  onto  eachother}  are  formed. 

The  potential  IFS-Codes  rn;  are  the  generating  transformations  of  the  largest 
orbits1*  of  each  class. 

Then  the  final  IFS-Codes  ©;  are  the  generating  transformations  of  the  largest  orbits 
within  each  C;  which  can  be  mapped  onto  eachother  by  (have  an  inter  orbit 
transformation  in)  {©j}+. 

Thus  the  final  and  also  minimal  IFS  can  be  calculated  by  the  following  iterative  process 

(a)  which  restrains  the  classes  to  those  orbits  which  can  be  mapped  onto  each  other 
only  by  m;  detected  so  far 

=>  more  and  more  classes  are  generated  (and  the  set  {m,}  is  increased) 

(b)  which  enlarges  the  classes  to  those  orbits  which  can  now  be  mapped  onto  each 
other  only  by  ©j  detected  so  far 

fewer  and  fewer  classes  are  generated  (and  the  set  {tn;}  is  decreased) 

This  process  is  repeated  until  no  more  (in  (a))  resp.  no  fewer  (in  (b))  classes  are 
generated.  Then  the  generating  transformations  of  the  largest  orbits  of  each  remaining 
class  form  a minimal  set  of  the  IFS-Codes  of  the  object. 

Since  a a-Hull  is  only  useful  to  find  the  orbits  on  the  boundary  of  an  object,  orbits  of 
inner  subobjects  (having  no  boundary  in  common  with  the  object)  cannot  be  found  in  the 
first  step.  But  if  you  form  the  difference  between  the  object  and  the  point  set  gener- 
ated by  the  IFS-Codes  in  the  first  step  one  will  get  a difference  object  Aj  in  the  second 
step.  Those  affine  transformations  mapping  boundary  orbits  (detected  by  different  a- 
Hulls)  of  the  object  A0  to  those  of  the  difference  object  A1  are  the  IFS-Codes  for  the  inner 
subobjects  now  lying  outside.  The  second  step  has  to  be  repeated  for  the  object  A,  (i  e N) 
as  long  as  the  difference  object  Ai+1  is  not  empty.  Adjusting  the  obtained  results  and 
eliminating  possible  deviations  by  comparing  the  inverse  IFS-mapping  of  the  hull  of 
inner  subobjects  with  the  object  hull,  also  the  inverse  IFS-Problem  for  overlapping 
subobjects  can  be  solved. 


* An  orbit  O is  a point  sequence  on  = 4>rrg(x)  VneN  where  $<n)  = f[  <k  <i>  affine  and  x e R2. 

n 

The  point  sequence  (orbit)  4>fni(p)  converges  for  a contractive  affine  maps  to  the  fixed  point  f of  <j> 
c { II  ®n  ffl,  e W u {e}  and  W+=  W*  \ {s}  with  W ={tBj  | cr  affine}  and  s the  Identity  Map. 

neN 

“ The  largest  orbits  of  a subobject  p(A)  (p  e W u {s})  are  those  orbits  O = {o0,  ...,on}  with  a 
minimal  contraction,  where  <j/n)  (o0)  e A for  n > 0 and  $(n)  (o0)  g A for  n < 0 and  <j>  e W+. 
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2 Definition  and  Application  of  the  O-Hull 

2. 1 The  a-HuIl  (fractal  hull)  in  Rx  <p 


Figure  1.  The  convex  Hull,  the  a-  and  cr-Hull  of  an  Object 


We  first  introduce  the  notion  of  generalized  spiral  disks.  Let  D*(f,  r)  denote  a spiral  disk 

00 

defined  by  the  point  set  [Jp,  with  fixed  point  (centre)  f,  the  non  negative  radius  r = Rekip 

1 

where  p;  e R x cp  and  d(f,  p,)  < r and  0 < <p  < 2mic,  m e N. 

A generalized  spiral  disk 


D’(f,  r)  = 


jD*(f,-r)c 

I D*(f,r) 


r <0 
r > 0 


(1) 


where  Xc  is  the  complement  of  X. 

Since  D'(f,  r)  is  not  solely  bounded  by  a log  spiral  (but  also  by  a part  of  the  q>  = 0 axis) 
a generalized  spiral  disc  segment  has  to  be  constructed  bounded  only  by  equally  curved 
log  spirals  having  the  same  fixed  point. 

Thus  a generalized  spiral  disk  segment  is  given  by  (Fig.  If. 


D(f,  r)  = D'(f,  Rekl?)  n D'(f,  -Re^71'^))  n D'(f,  Re^'L'P))  (2) 


Note:  For  negative  k one  will  get  a spiral  disc  segment  reflected  by  the  cp  = 0 axis  ! 

C0(S)  is  the  set  of fixed  points  of  spiral  disc  segments  of  radius  1 /a  that  have  S as  subset 

Co(S)  = {x€Rxcp|Sc  D(x,  ’/a) } (3) 


where  X is  the  set  closure  of  X. 

Let  X be  an  arbitrary  (finite  or  connected  set)  in  R x cp  then  the  intersection  of  all 
generalized  spiral  disc  segments  with  varying  fixed  point  x e X and  a fixed  radius  !/c 
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f Rx<p  X=0 

is  denoted  by  MC(X)  & n^.  1 / <t)  X*0  (4) 

[xeX 


where  0 is  the  empty  set. 

The  a-Hull  of  S is  the  intersection  of  all  generalized  spiral  disk  segments  of  radius  ]/a 
that  contain  all  points  of  S.  Thus  the  a-Hull  (fractal  hull)  is  given  by 


»a(S)  ^ Mo(Ca(S))  (5) 

A point  s e S is  termed  a-Extreme  (corner  of  »^o(S))  in  S,  if  there  exists  a D(c,  Vc) 
such  that  s lies  on  its  boundary  and  contains  all  other  points  of  S (see  Fig.  1). 

Thus  the  set  of  a-Extremes  is  defined  by 

Eo(S)  s {s  e S | 3c  € Cc(S):  s e 5(c,i/0)}  (6) 


where  oX  is  the  boundaiy  of  S. 

As  in  the  hmit  case  of  infinite  radius  (i.e.  Vk  close  to  0)  resp.  of  radius  = R (k  = 0)  the 
generalized  spiral  disc  segment  becomes  a half  plane  resp.  a circle  of  radius  R,  this  defi- 
nition includes  the  definition  of  the  convex  Hull  ^(S)  resp.  of  the  a-Hull  *?*a(S)  as  well. 


2.2  Defining  the  a-Hull  in  R2 


To  represent  a log  spiral  r = Rek<p  defined  in  the  polar  plane  (R  x <p)  an  equivalent  for  the 
x/y-plane  (R  x R)  has  to  be  found.  The  best  suited  equivalent  seems  to  be  one  formulated 
in  terms  of  affine  transformations.  In  addition,  this  representation  will  relate  log  spirals 
and  orbits  (lying  on  the  log  spirals)  which  union  forms  an  IFS  and  where  the  generating 
transformations  of  the  largest  of  these  orbits  are  the  IFS-Codes. 

Theorem  1 - Log  spiral  equivalent  - A log  spiral  r = Rekq>  can  be  represented  by  a 
continuos  function  of  powers  of  an  affine  transformation  m. 


Re**  = oj»  = R0(a)S(cxn,cyn)R0(-a)R0(n(-P)), 
where  n e R,  a e [-n,  k),  p e [-x,  n)  and  p * 0 


Lemma  1 Each  affine  transformation  can  be  represented  by  a symmetric  matrix  and  a 
rotation  matrix. 


Vra  = 


a b 
c d 


e R2x2  3p  e [-7C,  n)  3\y  = 


V 


(0  = V|'Rq(-P)  = 


r s 
s t 

( r cos/i  + s sin/3  s cos/3  - r sin  fr' 
scos/3  + t sin/3  t cos/3-s  sin/3 


j.  = aJ+tr+ad-bc 
■^(a  + d)2  +(c-b)2 

Proof  of  Lemma  1.  in 3 


t = 


^(a  + d)2  +(c-  b)2 


e R2x2  : 

with  tanp 
and  s = 


b-c 
a + d 

ac-bd 


^/(a  + d)2  +(c-b)2 


Lemma  2 Each  symmetric  matrix  can  be  represented  by  a rotation  a,  a scaling  matrix 
and  a rotation  -a. 
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Vvj/ : 


( r s^ 
Vs  0 


e R2x2  3 a e [-tc,  tc)  3 cx,  c e R : vj/  = R^a)  S(cx,  c ) RqC-cx) 


with  the  eigenvalues  cx 


r+t  + 


V(r-t)2 


+ 4s 


cy  = 


r+t-^(r-t)2  + 4s2 
2 


and  a = arctan-*—  the  angle  of  the  eigenvector  ex 


'c  -r> 


v s y 


Proof  of  Lemma  2.  in3 

Lemma  3 Each  affine  transformation®  maps  a point  x e Rek(p again  to  a point  x,  e Rek<p. 

Vm  e R2x2  and  x e Rek<p : rax  = x,  and  Xj  e Rek<p 
Proof  of  Lemma  3: 

r = Rek<p  resp.  ^/x2  +y2  = Re 


karctan- 


x = r cosep  = Rek<pcos<p  and  y = r sincp  = Rek'psin(p 


rax  = (» 


fxl 

| = 0) 

^Rekpcos  cp' 

k / acosf +bsin<jT 

Jy 

y Rc^  simp  j 

(ccoscp+dsin^ 

q>  = arctan--  = arctan  acos-?+  b sin  9 and 

X c cos  <p  + d sin  <p 


r = fx2  + y2  = Re^-y/Cacos  (p + b sin  (p)2  + (c  cos  (p + d sin  (p)2 


karctan 


c cos  <p+d  simp 


i — Kiuuaii — — i 

Rek<p>/(acos9+bsin9)2 +(ccos^+dsin9)  = Rek<(>e  acos^bsm^  j ;ek(p 

ccos<p+dsinp 


^/R2(acos<j£>+bsin<p)2  +R2(ccos<p+dsin<p)2  = Re 


k arctan- 


acosq)+bsin<p 

“V 4 


i o Ry  . y 

x-2  yz  arctan—^  = arctan— 

Rx  x 

Lemma  4 A power  n e G of  an  affine  transformation  maps  x e Rek,p  again  to  Xj  e Rek(f>. 

Vran  e R2x2  and  x e Rek<p : ©nx  = Xj  and  Xj  e Reklp  where  n e G 
Proof  of  Lemma  4: 

Lemma  3 holds  for  each  point  on  the  log  spiral  Rekcp 

I.e.  also  for  rax  and  for  ra2x  and  so  on  ...  until  ranx  resp.  also  for  ra_1,ra"2,  ...,ra'n 
Corollary  1 Each  power  of  an  affine  transformation  can  be  represented  as  follows. 


( a bY 


= Ro(«)S(cxn,  c/iRfll-alRoffil-P)) 

^c"cos2 a+c"sin2 a (c" -c")cosasina 
(c"-c")cosasina  c" cos2  a + c" sin2  a 


Ro(n(-P» 
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Proof  of  Corollary  1: 

Because  of  Lemma  4,  Lemma  1 can  be  also  defined  for  n e G : 

(t  s' 

con  = vj.cnR0(-p)n  where  vj/  = 

Vs  V 

^CD"=  (u'Rof-P))”  = rwr  = vfnRo(n(-P))  = 

Ro(«)  S(cx,  cy)  RflC-aiR^fa)  S(cx,  cy)  R^-cOR^a) ...  S(cx,  cy)  R^-oOR^nC-P))  = 
Rfl(a)  S(cx,  cy)  R^O)  S(cx,  cy)  R()(0) ...  S(cx,  cy)  Rfl(^)R<)(n(-p))  = 

Ro(a)  S(cx,  cy)n  Rof-^Rofnfa+p))  = R^a)  S(cx",  cy")  Ro(-a)  R^nf-p)) 

Corollary  2 Lemma  4 can  also  be  defined  for  powers  n e R instead  of  n e G. 

Vcon  € R2x2  and  x e Rek<f’ : conx  = Xj  and  x,  e Rek'fl  where  neR 
Proof  of  Corollary  2\ 

Lemma  4 also  holds  for  n e R,  if  con  is  represented  according  to  Corollary  1,  since  powers 
n € R not  applicable  to  <o  are  now  defined  for  cx  and  cy  (i.e.  n is  no  more  restricted  to  G) 

Summary  (Proof  of  Theorem  1) 

Lemma  4 - using  Lemma  3 - shows  that  ®n  (n  e G)  generates  point  sequences  on  log 
spirals  (the  so-called  orbits).  Corollary  2 points  out  that  for  n e R the  orbits  of  con  form 
log  spirals.  Corollary  1 defines  - by  means  of  Lemma  1 and  Lemma  2 - a representation 
of  ®n  = Ro(a)  S(cxn,  cyn)  Ro(-a)  R^nf-P))  which  can  be  used  to  calculate  con  in  R x R for 
neR,  what  concludes  the  chain  of  the  proof  for  Theorem  1. 


Definition  of  the  o-Hull  in  R2 

co 

Let  D*(f,  ®n)  be  a spiral  disk  (J  p,  with  fixed  point  f,  bounded  by  wn  with  0 < n < 2m x/p, 

i 

m e N and  radius  |<on(s-f)|,  where  p,  e R x R and  d(f,  p;)  < |wn(s-f)|,  s e R x R. 


Then  D’(f,  wn)  s 


|D*(f,ca“n)c 
1 D*(f,mn) 


n<0 

n^O 


(7) 


denotes  a generalized  spiral  disk  in  R2. 

The  generalized  spiral  disc  segments  described  by  means  of  log  spirals  in  R x cp  as  in  (2) 
can  now  be  defined  using  a continuos  function  of  powers  of  affine  transformations  in  R2: 

D(f,  o ) = D,(f,m,,)oD'(f,R0(7r/4)crn)nD,(f,R0(72)s(-l,l)ft)n)  (8) 


where  ron  = Rf)(a)S(cxn,  cyn)R0(-a)R0(n(-P)  and  0 < cx  < 1,  0 < cy  < 1,  -n  < a < n,  -k  < p < x 
Note:  for  negative  a and  p one  will  get  a spiral  disk  segment  reflected  by  the  x-axis 

Which  - by  adapting  the  definition  (3)  of  the  set  of  fixed  points  CCT(S)  to  R x R - leads  to 
the  definition  of  the  intersection  of  generalized  spiral  disk  segments  Ma(X)  resp.  the  a- 
Hull  in  R x R instead  of  R x <p  in  (4)  resp.  (5) 

[ RxR  X - 0 

MaPQ  = j P)D(x, 1 /|<un(s-f)|)  x^0  resP-  ^°(S) - Ma(Ca(S))  (9) 

[xeX 
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This  leads  to  additional  special  cases  for  the  a-Hull  defined  in  R2  : 
cx  = cy  in  R0(a)S(cxn,  ^"^(-aiR^ni-P) 

=>  D(f,  (Dn)  is  bounded  by  2 straight  lines  forming  an  angle  of  K/4 
If  cx  * cy  and  a = 0 in  R^a)  S(cx",  cyn)  R^-a)  R^nf-P)) 

=>  D(f,  ran)  is  bounded  by  3 exponential  parabolas. 

Both  shapes  are  necessary  to  find  all  orbits  (also  those  laying  on  non  spiral  like  curves) 
by  which  an  IFS  is  bound. 

2.3  Calculating  Orbits  using  the  a-Hull  in  R2 
Calculation  of  Orbit  Generating  Transformations 

Lemma  5 The  curvature  proportion  p of  ©n  = Rekc?  is  defined  by  a k independent  of  p 
Vos11  = Rek<p  with  p = R/r  : k = - 

2arctan 

Proof  of  Lemma  5 (cf.  Proof  of  Lemma  3) 

Using  v|/)T  = 


ln(c2  cos2  a+CySin2  a) 
(cx-cy)sin  acosa 
cxcos2  a+cysin2  a 


f r sN 

fXl 

Is  K 

V 

cx cos2  a+c"sin2  a (c"  -c")cosasina 
(cx  -c")  cos  a sin  a c"  cos2  a + c”  sin2  a 


for  the  substitution  of  x and  y resp.  r and  <p  in  the  definition  of  log  spirals: 
R^/lacoscp+bsin^)2  + (ccos<p+dsin<p)2  = Re 

n(cx-C 

c_  cos2  a+  c v sin2  a c ~ 

y for  cp  = 0 


X ( x^ 

of  Corollary  1 

UJ 


k arctan  C C0S  ^ ^ ^ 

acoscp+bsincp 


leads  to 


(cx+cycos2a)2  + ((cx— cy)  sina  cosa)2  = e 


. (cx-ctf) sin  acosa 
2karctan— — i— 


. (cx  cy)sinacos  i /■  *?  ->  « . « . 

2karctan — = ln(c2  cos2  a+c2  sin2  a) 


cxcos2  a + cvsin2  a 


Since  for  fixed  a (e.g.  tana  11)  different  values  of  cx  and  cy  form  log  spirals  Rekip  with 
arbitrary  k (different  in  p)  k resp.  1/k  can  be  represented  in  dependency  of  cx  and  cy: 


k = - 


2arctan 


^-y  rCSP-  k- 


2arctan 


(for  1/k  cx,cy  turn  into  and  !) 


Now  (A)  the  nominator  resp.  (B)  the  denominator  of  the  above  equations  for  k and  1/k 
can  be  equated  to  resp.  t0  calculate  cx  and  cy. 


2(k  + l) 


2(k  + l) 

-k?r 

(A)  =>  c2x+c2y  = 2e2Ck+l)  ..  ..  _x  _y  __  

(B)  =>  cx-cy  = -(cx+cy)tan^L_  for  k and 
This  leads  to: 


for  k and  c +cv  = 2e2(k  +1)  for  k"1 
x y 


for  k'1 
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q2  _ g2(k+l) 


k n 1 — tan 


4(k  + 1) 


~krc  1 + tan- 


1 + tan 


-k*1*  j 1-tan 
J^  = e2(k~'+l) 


4(k  + 1) 
n 


and 


C2=e2(k+')V 


4(k  + 1) 


1 + tan 


V 


4(k  + l) 


1 + tan2 


4^k  +1^/  an(j  _ e2(k'+!)  ^ 

> rr  1 * y 


1 + tan 


4(k_l  + 1) 


n 

4(k~*  + 1) 


1 + tan2 


4(k“‘  + l) 


fork 


fork'1 


(10) 


Note:  Now  it  obvious  why  k and  1/k  is  treated  differently,  because  otherwise  at  least  one 
of  the  resulting  cx  or  cy  can  become  greater  than  one  in  some  cases. 

The  representation  of  an  orbit  con  by  a continuos  function  of  n (Theorem  1)  allows  to 
form  the  a-Hull  for  an  attractor  A paramertized  only  by  cx  and  cy  and  a fixed  a (because 
in  that  case  an  arbitrary  p leads  to  equal  curved  c-Hulls  according  to  Lemma  5).  That  is, 
to  find  an  orbit  generating  transformation  co,  where  a maximal  number  of  discrete  bound- 
ary points  of  A lie  on  the  orbit  co"  (for  n € G)  forming  the  edge  of  a a-Hull  (for  n e R). 

Now  generating  transformations  cos  of  orbits  can  be  calculated  using  the  equivalent 
representation  of  a a-Hull  in  R2  if  a a-Hull  is  found,  where  at  least  four  boundary  points 
of  A pl5  p2,  p3,  p4  lay  on  one  of  the  circumscribed  spiral  segments  S with  the  fixed  point 
f,  forming  the  edge  of  this  (best  fitting)  hull. 

Thus,  Vpj  i e {1,2, 3, 4):  p,  € m"  4 ps  e R0(7c/4)co_n  4 pj  e R0(7t/2)S(-l,l)®n  and  con 
forms  the  edge  of  the  a-Hull. 

The  cx,  cy,  a and  p - for  the  log  spiral  equivalent  (Theorem  1)  - can  be  calculated  from 
cos  by  using  the  representation  eo-yR^-P)  and  v^=R0(a)S(cx,  cy)R0(-a)  (Lemma  1 and  2). 


Optimizing  the  orbit  generating  transformations  cos 

as  = R0(a)S(cxn,  cy")R0(-a)R0(n(-P)  can  now  be  optimized  in  the  following  way  : 

Since  p can  also  be  calculated  as  pk=Z(pkfR0(a(c)  Sfc'^c^y)  R0(-ak)(pk+]))  and  k=l,2,3 
=>  a new  p = (Pj+P2+P3)/3  can  be  obtained 
Now  a new  cx,  cy  and  a is  recalculated  using  pk  to  compute  yk: 


=>  cx  = (Cix+c2x+c3x)/3,  cy  = (cly+c2y+c3y)/3  and  a = (otj+a2+a3)/3 
rk  + tk  +ij(rk- 


where  ckx= 
for  v|/k  = oisRo(-Pk)  = 


■tk)2  +4sk 


V 


r^tk-3/(rkrtk)2+44  and  ^ =afctan 


(\  O 


U t 


k ) 


2 - sk 

: Ro(«k)  Cj.y)  R()(-ak)  and  k as  in  Lemma  2 


Reflection  transformations  (IFS-Codes)  - Either  cx  or  cy  is  negative 
If  two  best  fitting  a-Hulls  are  found  which  are  formed  by  the  log  spirals 
®,2"  (Pl)  = R^a)  S(cx2",cy2")  Ro(kx)  Ro(2n(-P))  or 
®22n+1  (P’i)  = Ro(«)  S(cx2n+1,cy2n+1)  Ro(-a)  Ro(-p(2n+l)) 
for  n e N°  resp.  1/n  e N and  p,,  p',  eR2,  one  affine  reflective  map  can  be  derived  : 
<Oi2=Ro(a)S(-cx,Cy)R0(-a)R0(-P)  resp.  w12=R0(a)S(cx,-cy)R0(-a)R0(-p), 
where  cx<0  resp.  cy<0,  if  R0(Y)o)1n(p,)=co2n(p'1)  resp.  R0(-y)co1n(p1)=(B2n(p',)  for  y>0. 
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2. 4 How  to  calculate  IFS-Codes 
Deriving  the  IFS-Codes  from  boundary  orbits 

Form  larger  and  larger  curved  cr-Hulls,  where  the  cx  and  cy  are  computed  according  to 
Lemma  5 by  using  smaller  and  smaller  k resp.  1/k  (for  Rek'f)  in  the  log  spiral  equivalent 
of  Theorem  1)  for  a fixed  a (Fig.  2). 

Each  time  you  get  a best  fitting  hull  (at  least  four  boundary  points  of  an  attractor  A lie 
on  the  edge  of  the  cr-Hull)  the  ®s  - generating  that  orbit  defined  by  these  4 points  - is  one 
IFS-Code  of  the  IFS  Wm  for  a subobject  lying  outside  (i.e.  having  a part  of  its  boundary 
in  common  with  the  object)  as  long  as  ®s  4 ®s,  (Fig . 2).  I.e.  Wi+1  = W;  u (®s)  for  i = 
{1,  . . .,  m},  if  ®s  4 ®s.  where  ®s,  e Wj  and  cos  = ®s,,  if  ®s,  = rosros»  where  cos*  e W; 


Figure  2.  IFS-Code  calculation  by  orbit  generating  transformations 

Besides,  this  extremely  simplifies  the  classification  process  (cf.  introduction)  because 
larger  orbits  are  found  first  and  the  smaller  ones  have  not  to  be  put  into  the  classes  any 
more.  Using  wn  = RfJ(a)S(cxn,  cyn)R0(-a)R0(n(-(I),  several  special  cases  can  be  apriori 
excluded  and  thus  unnecessary  classes  are  not  formed: 

• Multiples  of  other  orbits  are  detected  if  S(c’x,  c'y) = S(cxn,  cyn)  n e G resp.  if  as  = as.. 

• The  orbit  generating  transformations  cos  of  the  largest  orbits  of  one  class  is  that  w's 
with  the  largest  p. 

• Orbits  which  can  be  affinely  mapped  (by  inter  orbit  transformations)  onto  eachother 
will  now  be  detected  by  those  cos  having  a multiple  in  p - without  using  an  affine 
invariant  representation7-8. 

Thus  in  the  classification  process  it  is  not  longer  necessary  to  use  classes  of  orbits  but 
classes  of  the  orbit  generating  and  inter  orbit  transformations  represented  by  their  log 
spiral  equivalent  (of  Theorem  1) 

Deriving  the  IFS-Codes  from  non  boundary  orbits 

Find  all  rokj  that  for  at  least  one  of  the  spirals  Sj  s (which  are  formed  by  each  best  fitting  a 
-Hull  of  A for  the  IFS-Codes  ©Xi)  circumscribed  to  A\ro^'kj  (U®k-n  A))  holds  o)kj(Si;)  = 

i 

Ski  and  Ski  a circumscribed  spiral  to  A \ y<ok_H  (A)  and  (okj  the  IFS-Codes  for  an  inner 
subobject  (i.e.  having  to  less  boundary  points  in  common  with  the  object). 
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i 

Repeat  this  process  until  A\|J<uki(A)  = 0,  where  {<okj}  ={(ok.lj}u{cokj}  and  k ={2,...,n}. 
Note:  Thus  the  IFS-Codes  of  attractors  with  overlapping  subobjects  can  also  be  calculated. 


3 Conclusion 

3. 1 Summary 

In  this  work  is  shown  how  the  boundary  of  a discrete  self  affine  object  can  be  calculated. 
This  is  the  proposition  to  form  orbits  of  boundary  points  of  an  object.  Subsequently  the 
orbits  have  to  be  classified  to  calculate  the  IFS-Codes  (affine  maps)  of  a minimal  IFS  as 
generating  transformations  of  the  largest  of  these  orbits  in  each  class. 

The  boundary  orbits  are  found  by  using  a generalized  convex  hull  - the  cr-Hull.  This 
hull  is  formed  by  the  intersection  of  generalized  spiral  disk  segments  (parts  of  the  plane 
bounded  by  3 log  spirals  equal  in  curvature  proportions  and  fixed  points)  instead  of  half 
planes  used  for  the  convex  hull.  To  reduce  the  calculation  costs,  an  equivalent 
representation  of  spirals  (in  the  polar  plane)  is  defined  by  powers  of  affine 
transformations  (in  the  x/y-plane). 

This  representation  has  the  additional  advantage,  that  - by  the  used  affine 
transformations  - orbits  can  easily  be  calculated  for  the  subsequent  classification. 
Besides,  the  affine  equivalent  to  log  spirals  extremely  reduces  the  expense  for 
comparison  and  differentiation  of  the  orbits  within  the  classification  process. 

Further  more  generating  transformations  for  non  boundary  orbits  can  be  detected  after 
removing  the  subobjects  lying  outside  which  have  been  already  calculated.  Thus  the  IFS- 
Codes  for  subobjects  having  no  boundary  in  common  with  the  object  can  be  computed, 
even  if  the  subobjects  are  overlapping. 

3.2  Future  Work 
Implementation 

This  is  a theoretical  work,  which  forms  the  basis  for  the  IFS-Code  calculation  and 
classification  resp.  overcomes  the  difficulties  of  their  calculation  existing  so  far. 
Therefore  an  implementation  will  be  made  together  with  the  elaboration  of  the  classi- 
fication part  which  has  to  be  completed  to  solve  the  inverse  IFS-Problem  as  a whole. 

But  the  prototype  developed  for  the  last  work9  already  shows,  that  in  principle  this 
problem  is  solvable  by  using  a less  general  hull  similar  to  the  cr-Hull  and  by  mapping 
orbits  of  boundary  points  onto  eachother. 

Besides,  the  experience  made  by  the  implementation  of  this  prototype  showed  the 
need  to  have  a more  general  tool  (the  a-Hull)  to  calculate  and  classify  the  orbits  as  basis 
for  the  computation  of  a minimal  IFS. 

Analysing  IFS-attractors  under  the  aspects  of  the  O-Hull 

So  far  the  affine  equivalent  of  log  spirals  (the  basis  of  the  a-Hull)  has  primarily  be  used 
to  calculate  the  boundary  of  an  IFS-object.  But  it  can  also  be  an  instrument  to  categorize 
EFS-objects  according  to  the  intrinsic  properties,  such  as  structure,  location  and  form  of 
the  object  parts  (subobjects)  resp.  orbits. 

Therefore  the  correlation  between  orbits  and  the  parameters  of  the  a-Hull  such  as 
scaling  factors  (and  powers  of  them)  - determining  the  curvature  of  the  orbits  - and  the 
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two  angles  - one  determining  the  distance  of  orbit  points  and  the  other  determining  the 
contraction  of  an  orbit  - have  to  be  studied  in  further  detail. 

The  gained  knowledge  will  also  improve  the  IFS-calculation  and  contribute  to  a better 
adaptation  of  the  a-Hull. 

Can  the  classification  process  be  more  improved  by  the  a-Hull  ? 

Some  implications  of  the  parameters  of  the  a-Hull  for  the  classification  have  been 
already  analysed  in  this  paper,  but  there  should  be  more  implicit  properties  of  the  a-Hull 
which  can  be  used  to  classify  the  orbits. 

Besides  the  choice  of  the  order  - not  only  in  dependency  of  a decreasing  function  of 
the  curvature  proportion  - by  which  the  a-Hulls  are  formed  can  influence  the 
discrimination  process  of  the  detected  orbits  to  be  classified. 

Last  but  not  least  the  affine  equivalent  to  log  spirals  can  be  used  - within  the 
classification  - to  decide  if  one  matrix  can  be  divided  by  another  resp.  for  the  comparison 
of  matrices  with  respect  to  their  fixed  points. 

Designing  a more  adequate  a-Hull 

Based  on  the  present  knowledge  about  orbits  of  boundary  points  of  IFS-objects,  the  basic 
element  of  the  a-Hull  was  designed  as  the  intersection  of  3 spiral  disks  equal  in 
curvature  proportion  and  fixed  point.  But  in  praxis  it  may  turn  out  that  more  or  less  log 
spirals  with  different  curvature  and  fixed  point  will  do  a better  job. 

This  variations  also  will  change  the  size  of  the  basic  element,  what  will  influence  the 
number  and  accuracy  of  orbits  found  on  the  boundary.  That  is,  roughly  speaking  how 
"tight"  a-Hulls  can  be  formed. 

So  far  only  two  types  of  basic  elements  have  been  investigated.  An  intersection  of  two 
log  spiral  disks  having  the  same  curvature  and  the  same  fixed  point.  But  - though 
simpler  to  handle  - this  leads  to  basic  elements  infinite  in  size,  which  are  not  sufficient  to 
circumscribe  an  IFS-object  in  several  cases. 

Another  basic  element  was  formed  by  two  log  spiral  disks  identically  curved  in 
opposite  direction  and  having  the  same  tangent  in  their  two  tangent  points.  This  leads  to 
more  consistent  (with  respect  to  the  definition  of  hulls)  mathematical  model  of  the  a- 
Hull,  but  this  basic  element  is  convex  and  using  it's  complement  again  leads  to  less 
"tight"  a-Hulls,  which  prevent  that  the  fixed  points  of  orbits  to  lay  on  the  boundary. 

Generalization  of  the  a-Hull 

In  general  each  discrete  object  can  be  represented  by  an  IFS  as  long  as  sufficient  many 
IFS-Codes  are  used  (in  the  extreme  case  one  IFS-Code  per  pixel  !).  Decomposing  an 
object  in  adequate  parts  will  reduce  the  number  of  IFS-Codes10.  So  the  question  is  not  to 
design  more  and  more  complex  hulls  (bounded  by  highly  non  linear  curves11),  but  to  use 
the  a-Hull  to  find  an  minimal  number  of  orbits  fitting  best  into  the  calculated  boundary 
points.  Two  possible  solutions  are  outlined  below: 

One  to  enlarge  the  pixel  size  by  an  appropriate  scale  and  to  fit  the  orbits  into  less 
boundary  points.  Another  to  make  the  pixels  smaller  - so  that  each  point  becomes  a 
boundary  point  and  to  find  all  different  orbits  within  this  unconnected  point  set.  The 
second  approach  seems  to  be  more  complicated  (complexity  in  the  point  relation),  but  it 
solves  the  encoding  problem  in  one  step,  because  there  are  no  subobjects  lying  inside, 
anymore. 
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